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I. INTRODUCTION 

Shell model is a recognized tool for microscopic stud- 
ies of nuclear structure. No-core shell model (NCSM) 
Q 13 1 the version of the shell model where all A nu- 
cleons are spectroscopicallr active, is widely used now 
(see, e. g., Refs. 0)11: IS 111) for ab initio calculations of 
light nuclei (up through A = 16) with modern realistic 
nucleon-nucleon and three-nucleon forces. 

NCSM utilizes the basis of Slater determinants of 
single-particle oscillator states. These basis functions are 
known to have spurious contributions of center-of-mass 
(CM) excitations. The wave functions of physically ac- 
ceptable eigenstatcs of intrinsic NCSM Hamiltonian 

1 ^ (p_p.)2 A A 

Ha^ + XI ^W^f.y + H VArjvAr^yfe, 

i<3 i<j i<j<k 

(1) 

where m is the nucleon mass, Vmnaj is the two-nucleon 
interaction (including both strong and electromagnetic 
components), VNNN,ijk is the three-nucleon interaction, 
should be arranged as spurious-free linear combinations 
of basis states. 

To achieve this, the auxiliary Hamiltonian 

Hncsm = Ha + PQo (2) 

is conventionally diagonalized within NCSM instead of 
the Hamiltonian (Q. Here 

go = HcM - \h^, (3) 

HcM — TcM + UcM (4) 

is the harmonic oscillator CM Hamiltonian, Tcm is the 
CM kinetic energy operator, and 

UcM = ^^Amn'B?, (5) 

where 

R=lf]r,. (6) 

The term /SQo with large enough parameter (3 has no 
effect on the intrinsic states of the ^-body system due 
to the translational invariance of the Hamiltonian 
it shifts up in energy spurious CM-excited states and 
projects out the spurious contributions in the low-lying 



eigenstatcs. As a result, the physical low-lying eigen- 
statcs of (121) correspond to the QMl CM-excitation and 
are independent on the choice of (3. 

I suggest below a projection operator Pcm that can be 
used to project out spurious CM-excited components and 
to obtain spurious- free linear combinations of basis Slater 
determinants that can be used as a new spurious-free ba- 
sis for direct diagonalization of the intrinsic Hamiltonian 
(|H). The complete spurious-free basis corresponding to 
the QhVt CM-excitations, is much smaller than the basis of 
all Slater determinants including all kHVL CM-excitations 
with K < N where N is the maximal oscillator quanta 
of the NhVl NCSM model space used in the calculations. 
Therefore it is expected that the use of the projection 
operator Pcm will simplify essentially the NCSM stud- 
ies of nuclear structure, will make it possible to arrange 
the calculations in a larger NhVl model spaces with the 
same computer facilities and hence to improve the accu- 
racy of the NCSM predictions, etc. 

I note also that the so-called m-scheme is convention- 
ally utilized in the NCSM applications, i. e. the basis 
Slater determinants are used that do not have definite 
values of the orbital angular momentum L, of the total 
angular momentum J, and of the total spin S. The m- 
scheme makes it possible to use well-developed computa- 
tional methods and available respective computer codes. 
However the basis of the m-scheme Slater determinants 
is very large since it includes all states with all possible 
values of J < Jmax, L < Lmax and S < S'max where the 
maximal values Jmax 7 -^max and S'max are large enough in 
modern NCSM applications and depend on the particu- 
lar nucleus under consideration and on the Nhfl model 
space used in the calculations. 

I suggest below the projection operators Pj, Pl and 
Ps on the states with definite J, L and S values. These 
projection operators as well as Pcm can be easily utilized 
within the existing NCSM codes to reduce essentially the 
number of the basis states. 



II. CM-PROJECTOR Pcm 

Let 

N 

* = ^ a^^^ (7) 

be a vector (wave function) defined in the Nhfl model 
space; N is the maximal possible CM-excitation quanta 
in this model space. Equation Q presents expansion 



2 



of ^ in the series of functions with a definite CM- 
excitation quanta x = 0, 1, N . The functions ^ are 
the eigenfunctions of the harmonic oscillator CM Hamil- 
tonian: 



(8) 

Due to Eq. ©, the operator Qo defined by Eq. © 
acts as 'anti-projector': it projects out the spurious- free 
component ^'s/ = ^'o of the wave function, 

Qo'i'sf = 0. (9) 
We can also define anti-projectors 

Q^ = HcM - (^x+l^nn (10) 

which project out components with given values of the 
CM excitation quanta k: 

Q^^^^Q. (11) 

We can extract the spurious-free content of '5 by 
the subsequent use of the operators ifTnji : 



*2 = Q2*l, 



(12a) 
(12b) 



^sf = *w = Qjv*jv-i. (12c) 
Equations p2|l are equivalent to the following equation: 



where the operator 



N 



(13) 



(14) 



Let us call P 'quasi-projector'. Mathematically P is 
not a projection operator since it does not fit the stan- 
dard property of the projection operators, 



= P. 



(15) 



The function is a superposition Q of the spurious-free 
sj = ^0 a-nd spurious components with >r ^ 0. The 
standard projection operator property (|15|l guaranteers 
that 



(16) 



Instead of (|16|l . the quasi-projector P when applied to ^' 
results in 



(17) 



The constant D can be easily calculated using Eqs. 
and iCnil: 



D 



N 



(18) 



To become a projector, the quasi-projector P should be 
'normalized': 



(19) 



In applications, one can use either P or P. Really, it is 
usually needed to extract from ^ its normalized spurious- 
free component ^'s/- The multiplier ao is usually un- 
known. Hence after using either the quasi-projector H14|l 
or the projector H19|) . one needs to normalize either the 
function Dao'^sf or the function ao^'s/. Clearly, the 
same computational efforts are required to normalize the 
functions Dao'I's/ and ao^'s/- 



III. OTHER USEFUL PROJECTORS 

The same idea can be utilized for the construction of 
other useful projectors. As an example, let us construct 
the projector on the states with a definite value of the 
angular momentum. 

Let = L'^ + Ly + be the standard orbital momen- 
tum operator. Its eigenvalues are known to be L{L + 1). 
We define now the operators 



}L = L^~LiL+l) 



and 



L-l 



P^ 



uq^ n 



(20) 



(21) 



K = L+1 



tum in the given Nhfl shell model space. The non- 
normalized quasi-projector H21|l can be used like the 
CM non-normalized quasi-projector H14II to extract (non- 
normalized) component with the definite value of the or- 
bital momentum L by the algorithm described briefly by 
Eq. (221 or in more detail by Eqs. ifT^ . 
The projector p^™'™ can be expressed as 



P^ 



1 



-P^ 



(22) 



where 



L-l 



DI--^\{[k{k + 1)~L{L + 1)] 



K=0 



X n + + (23) 
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The structure of the projectors pj^'^^ ^ pS^n^ ^-^^ 
pTmax Qj-^ ^Yie states with given values of the total an- 
gular momentum J, total spin S or isospin T, is exactly 
the same. The only difference is that in the case of an 
odd-A system, one should use half-integer J, T oi S val- 
ues and modify respectively the products in Eqs. pifl 
and if^. 

The standard projection operator property H15|l is valid 
for all projectors (but not quasi-projectors) discussed 
above. 



IV. CONCLUSIONS 

Expression of the projection operators on the states 
with definite value of the angular momentum in the form 
of the expansion in the powers of the SU(2) generators, 
are known in the literature (see, e. g., Ref. Q)- However, 
in the general case, this polynomial includes an infinite 
number of terms and is inconvenient for the use in the 
nuclear shell model applications. As it was shown above, 
in the case of the shell model, the projector can be taken 
in the form of a finite polynomial in generators that is 
much more useful for the applications. The suggested 
projectors p^'^'^'^^ pJmax^ pSm^x g^jj^ p^'^'^^ are of this 
form. 

The CM-projector Pcm is also suggested as a finite ex- 
pansion in the powers of a simple CM harmonic oscillator 
operator Hcm- To the best of my knowledge, similar ex- 
pression for the CM-projector were never discussed in the 
literature. 



The Lanczos iteration approach is utilized in the mod- 
ern shell model codes, i. e. the basis vectors are obtained 
successively by acting by the Hamiltonian on the vector 
obtained on the previous step. The intrinsic Hamiltonian 
(UJ and the NCSM Hamiltonian (0) cannot produce CM- 
excited states or to change the value of the total angular 
momentum of the state. Hence it is possible to project 
only the pivot vector (the initial vector in the Lanczos it- 
eration approach) on the spurious-free subspace with the 
given definite value of the total angular momentum J; all 
the rest basis vectors will be produced spurious-free and 
with the same value of J by the Lanczos iterations. 

Formally one can use the projected pivot vector and 
the intrinsic Hamiltonian Q instead of the auxiliary 
Hamiltonian ^ in the NCSM applications. However it 
is well known that the spurious states will be produced 
in the Lanczos iteration approach due to the computer 
noise (round-off errors). The term PQq in Eq. (O stabi- 
lizes the NCSM calculations reducing the computer noise 
if /3 is sufficiently large. Therefore it looks reasonable to 
utilize the auxiliary Hamiltonian in the applications; 
probably it is reasonable to add the term 7I — J( J-|- 1) 
with sufficiently large 7 to the Hamiltonian to reduce 
the computer noise in the calculations of the states with 
the definite value J of the total angular momentum. 
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